An 77-condensate of fermionic atom pairs via adiabatic state preparation 
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We discuss how an ^-condensate, corresponding to an exact excited eigenstate of the Fermi- 
Hubbard model, can be produced with cold atoms in an optical lattice. Using time-dependent 
density matrix renormalisation group methods, we analyse a state preparation scheme beginning 
from a band insulator state in an optical superlattice. This state can act as an important test case, 
both for adiabatic preparation methods and the implementation of the many-body Hamiltonian, 
and measurements on the final state can be used to help detect associated errors. 
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Experiments with cold atoms in optical lattices not 
only make possible the realisation of many-body lattice 
Hamiltonians and their corresponding ground states [l], 
2], but also exhibit long coherence times. This opens the 
way to produce excited many-body states and consider 
the related quantum dynamics, as demonstrated by re- 
cent investigations of repulsively bound atom pairs 0, 0] • 
A key question in this context is how to prepare specific 
excited states, especially those corresponding to inter- 
esting quantum phases. Here we show that exact ex- 
cited eigenstates of the Fermi Hubbard model, the 77- 
condenstates first discussed by Yang [5| can be realised 
in experiments by combining an adiabatic ramp begin- 
ning from an insulating state in an optical superlattice 
with a sudden switch in the interaction strength (see 
Fig. la). These states exhibit long range order in all di- 
mensions and have been discussed in the context of high 
temperature superconductivity Q. Moreover, as exact 
excited eigenstates they provide (i) an ideal test case for 
the use of adiabatic ramping processes in state prepara- 
tion which has important possible applications in 
the production of low-entropy ground states, and (ii) the 
possibility to validate the implementation of the many- 
body Hamiltonian, by testing the properties of the final 
state. 

Below we show that the state preparation process pro- 
ceeds with high fidelities for realistic experimental size 
scales and parameters, even in the presence of imper- 
fections and noise. We focus on the ID case, where 
time-dependent density matrix renormalisation group 
(TDMRG) methods [10] allow exact calculations for rel- 
evant experimental conditions. However, the properties 
of the fy-condensate are essentially identical in higher 
dimensions, and we expect that this switch and ramp 
scheme will work similarly in 2D and 3D. We also show 
that the superlattice scheme has strong advantages over 
an alternative schemes involving the adiabatic opening of 
a harmonic trap Q . We then discuss how errors in state 
preparation or implementation of the Hubbard Hamilto- 
nian can be revealed and characterised in experiments 
via measurements made on the ry-condensate. 

The target state of our switch and ramp process, the 
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FIG. 1: a) Preparation of an ^-condensate: 1) Begin in an 
insulator state \ipi) with attractive onsite interactions U in an 
optical superlattice with depth Vsl', 2) Switch U to a positive 
value larger than the bandgap; 3) Delocalise onsite pairs by 
adiabatic removal of the superlattice. b) Full spectrum of 
energies for Hfh with 17 > in ID for a single pair of atoms, 
one of each spin species, plotted as a function of centre of 
mass quasimomentum k, with a the lattice spacing. 



77-condensate, is an exact excited eigenstate of the Fermi 
Hubbard Hamiltonian (h — 1) in D dimensions 



H 



FH 



u 



(1) 



This Hamiltonian describes the dynamics of atoms in 
the lowest band of an optical lattice [3, [ll|, with c\_ a 
a fermionic annihilation operator for particles of spin 
a {t, 1} on lattice site i = (ii . . . , id), J the tunnelling 
amplitude, U the onsite interaction energy shift, and 
m tCr = c\ a Ci a . The i] condensate can be constructed via 
the operator = J2i(~ 1)^° =1 td c\ i-Cj j nrs t introduced 
by Yang, which has the property [HpH,rf] = Urf . The 
state \rjjsr) ~ (rf) N |vac) is an eigenstate of Hfh with 
energy NU for positive integer N. Below we focus on 
the case U > 0, where \t]n) is a condensate of N re- 
pulsively bound atom pairs [3j. In Fig. lb, we plot the 
eigenenergies of Hpu when we have one particle of each 
spin on a ID lattice, as a function of the centre-of-mass 
quasimomentum. The single r\ pair is indicated in the 
plot, and corresponds to a repulsively bound onsite pair 
at the edge of the Brillouin zone, i.e., with center-of-mass 
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quasimomentum n/a. 

Switch and ramp process:- The 77-condensate with N 
pairs can be prepared using a switch and ramp process, 
combining an adiabatic ramp with a sudden switch in the 
interaction strength. Adiabatic ramps have previously 
been discussed for preparation of many-body ground 
states in optical lattices 0, S|- I n an adiabatic ramp, 
one prepares a state \ipf) of a Hamiltonian Hq beginning 
from a non-degenerate, gapped initial state \tpi) that is 
an eigenstate of the Hamiltonian Hq + V. By removing V 
adiabatically, the state follows the instantaneous eigen- 
states of H + V(t) and ends in \tpf). The key is that 
\t/ji) should be a gapped state of Hq + V that is easy to 
prepare with very low entropy via standard cooling and 
loading techniques @, EH ■ 

Here we propose to begin from a band insulator in 
the lowest sites of an optical superlattice [HI], as de- 
picted on the left in Fig. 2a, which is the ground state 
of Hfh + V, with V the Hamiltonian describing the su- 
perlattice potential. For the case depicted in Figs, la, 
2a, where the superlattice period is twice the original 
lattice spacing, V = VsL^2i even ni- This state has an 
energy gap €sl ~ Vsl corresponding to the superlattice 
bandgap, and a filling factor which is set by the super- 
lattice period Q (e.g., half filling in Figs, la, 2a). If wc 
were to let Vsl adiabatically we would connect this 
ground state to the ground state of Hfh- Instead, we 
can suddenly switch U (on a timescale short compared 
with J^ 1 ) to a value larger than €$l (see Fig. la). In 
the limit \U — €sl\ 3> J, this switching will create an ex- 
cited eigenstate of Hfh + V, as shown in the transition 
from the left panel to the right panel in Fig. 2a. Adi- 
abatic removal of the superlattice, Vsl — ► 0, will then 
lead to an excited eigenstate of Hpu- This latter state 
will correspond to the lowest energy state in which all 
particles exist in repulsively bound pairs, which is the 
77-condensatc. The energy spectrum of the Hamiltonian 
for a small system in ID is plotted during the ramp in 
Fig. 2a, and we see that the state is always separated by 
a gap ~ U to lower lying states, and the gap from the su- 
perlattice ~ Vsl. Whilst in general, the adiabatic ramp 
could be optimised using optimal control methods [lj], 



we choose a simple exponential ramp for the superlat- 
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the approximate linear dependence of the gap on Vsl- 
Here, T is the total ramp time, and v the ramp speed. 
Note that for large system sizes, the energy gap can be- 
come very small, but that for finite systems, a gap will 
always exist to the other excited states. The key question 
is how slow this ramp should be in order to obtain the 
^-condensate with high fidelity for realistic system sizes 
- 100 sites [l5 1 . 

Fidelity measures:- We measure closeness of the final 
state \ipf) to the 77-condensate in two ways: a) Via the 
full many-body fidelity T = \(ipf\ilN)\ 2 , and b) via the 
similarity of characteristic correlation functions of \tpf) 
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FIG. 2: a) Energy eigenstates of Hfh + V as a function of 
time for a small example system with iV| = Ni = 3, M = 
6. Left: Lowest energy states for strong initial attraction 
U/J — —30, in the presence of a superlattice. The lowest 
energy state is the initial state in our preparation scheme, 
\i/>i). The shaded area denoted the excited manifold of states 
with one dissociated pair. Right: The highest energy levels of 
Hfh + V(t) as a function of time during the adiabatic ramp. 
The lowest energy eigenstate \ip(t)} in the upper manifold 
where all atoms exist in pairs is equal to |t7jv) at t = T. During 
the ramp, a gap of order U always exists to the manifold 
(shaded area) where some atoms are unpaired, and a gap to 
higher levels is present when the superlattice is present, b) 
Pair momentum distribution of a perfect 77-condensate 

(see text). 



to those of the 77 state. Remarkably, we will show below 
that fidelities T ~ 1 can be obtained for long ramps, de- 
spite the fact that this quantity is exponentially sensitive 
to the system size, due to the increase in the size of the 
many-body Hilbert space. Indeed, we note that in large 
systems, states close to \t]n) can have essentially the same 
physical character as the desired state, and the associated 
correlation functions may not be significantly changed by 
a few small defects in the state, even if T becomes small. 
We thus also consider the comparison between charac- 
teristic correlation functions for the final state and ^at), 
which gives a measure that can be directly measured in 
experiments, and is not eponentially sensitive to the size 
of the system. In particular, we are interested in the pair 
momentum distribution C&(t) = Ci , which can be 
measured, e.g., by associating atoms in doubly occupied 
sites to molecules, and releasing them from the lattice to 
perform a time-of- flight measurement. This correlation 
function is strongly peaked for l?^}, reflecting the off- 
diagonal long-range order (ODLRO) exhibited by the 77- 
condensate in any dimension, with the pairing correlator 

Cj^fj = (w|cL T c LiC„ac„, T |77Ar) - hi, N e^ m - n y M , 
(if m ^ n), and I m ,n = N(M - N)/(M - 1). The pair 
momentum distribution is the Fourier transform of this 
quantity, C^ n) = £ r 
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(see Fig [5b). We will also consider the total distribution 
distance V{t) = l-£ k |C k (i)-4" N> |/ £ k |C k (i)+c[" N> |. 
Many-body Fidelities:- In Fig. 3a we plot the fidelity 
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FIG. 3: a) Fidelities for the superlattice and parabolic trap 
ramp as a function of ramp time T, computed using Hxy- 
The superlattice ramp shape is Vsl = 1J{e~ vt — e~ vT )/(l — 
e~ vT ), v = J/8. For the parabolic trap, we use the same 
shape with initial Vp/J = 0.1, v — J/12. U is decreased with 
the same shape as the potential in each case, with U = 30J 
at t — 0. The inset shows results for longer ramp times 
with M = 16. b) Onsite pair momentum distribution after 
T — 2400 J -1 for the superlattice (solid lines) and parabolic 
trap (dotted lines) ramps for M — 32 (left) and M = 64 
(right), computed using Hxy- c) Final state fidelity T as a 
function of correlation function distance T>(T) from the per- 
fect 77-condenstate, computed using Hxy- d) Cfc(T) for su- 
perlattice ramps, computed using Hfh, with a number of 
impurities Ni = 1 (left) and N = 2 (right), for T = 200J" 1 
(solid black), and T = 400J" 1 (grey dashed). 



T at the end of the ramp as a function of ramp time T 
for different system sizes M. In order to perform more 
accurate calculations for reasonable computational time, 
these results are obtained in the limit U ^> J. On states 
that have only repulsively bound pairs, Hamiltonian ([I]) 
acts as H XY = —(J 2 /U) S i S j + 2V SL £ ievcn Sf in 

second order perturbation theory, with S, = (Sf , Sf , Sf) 
denoting a vector of spin-1/2 operators, and spin states 
corresponding to sites that are occupied or unoccupied 
by a pair of atoms [l^ | . Remarkably, for long ramp times 
it is possible to obtain unit fidelity, i.e., essentially per- 
fect ^-condensates. The fidelities are also high for typical 
experimental sizes and shorter timescales, with M = 64, 
T < lOOOJ -1 . Although the timescales required to ob- 
tain a fixed fidelity increase with system size, we note 
(i) that we are already in the regime of experimentally 
relevant system sizes, and (ii) that the sensitivity of T 
increases exponentially with the size of the system, as 
discussed above. 

Pair momentum distributions:- This picture is comple- 



mented by the pair momentum distributions, depicted in 
Fig. 3b. In each case, the Ty-pairing peaks C„/ a are clearly 
visible, though for ramps with final fidelity lower than 
one, these peaks are somewhat broadened. In Fig. 3c we 
quantify this relationship between the fidelity T and the 
overlap of the pair momentum distribution with that of 
the perfect ^-condensate, as measured by T>{T). Over a 
wide range of T and for different system sizes, we see that 
these quantities are strongly correlated, so that sharpness 
of the peak could be used to infer the quality of the 77- 
condensate in experiments. 

Comparison with opening a harmonic trap:- For the 
same range of T and M = 32,48,64 we also com- 
pare our superlattice scheme to an adiabatic prepara- 
tion scheme that was recently proposed, in which a band 
insulator is formed in the centre of a harmonic trap 
Vtrap = J2i ^p(* a ) 2 j an d the trap is then opened to pro- 
duce the final state @. As shown in Fig. 3a, we see that 
for the same system sizes and ramping times, we obtain fi- 
delities that are roughly two orders of magnitude smaller 
from ramping the harmonic trap. For M > 32, we see 
poorer scaling for the harmonic trap ramps than for the 
superlattice ramp (For 64 lattice sites we obtain fideli- 
ties T ~ 10~ 12 ). Similar effects are seen in Fig. 3b in 
the broadening of the final pair momentum distribution. 
The superlattice ramp appears to have an advantage over 
the harmonic trap scheme because the atoms do not need 
to tunnel across the whole system during the ramp, but 
rather establish coherence locally. 

Imperfections:- We now investigate imperfections in 
the state preparation process. We will start by address- 
ing how missing atoms in the initial state, noise, and 
harmonic trapping potentials affect preparation of the 
?7-condensate. We will then discuss how measurement 
of time-dependence of correlation functions for the final 
state can be used to reveal and characterise these imper- 
fections in experiments. 

Imperfections - missing atoms:- To study the impact of 
missing atoms in the initial insulator state, we computed 
the time-evolution of the adiabatic ramp (with the full 
Hamiltonian) starting with localised defects. Regardless 
of where these defects are present, and whether we have 
only missing atoms or complete missing pairs, this re- 
sults in a broadening of the peaks in the pair momentum 
distribution. Examples are shown in Fig. 3d for a ramp 
at half filling with a number of missing atoms N% = 1,2. 
The resulting correlation functions are, however, stable 
in time (see below for further discussion). 

Imperfections - noise:- Motivated by recent discussions 
[ijj], we also investigated this ramp in the presence of 
noise. This would primarily arise from fluctuations in the 
lattice depth, which would change the value of J. Note 
that in the superlattice ramp, J (coupling neighbouring 
sites) is always non-zero, even though the effective tun- 
nelling at the beginning of the ramp is made small by 
the superlattice, ~ J 2 ' jtsh- With a variation of J up to 
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FIG. 4: Stability of a state close to an 77-condensate with 
imperfections, a) Time dependence of C v / a for the initial 
state with impurities as defined in the text on 32 sites, with 
16 onsite pairs and varying impurity count Ni, for U/J = 4 
(dashed black), U/J = 10 (solid grey), b) Same as (a), but 
with additional trapping potential Vp/J = 1.25 x 1CP 4 . 



10% with a variety of correlation times for the noise, we 
found no significant effect on the final state fidelity. 

Effect of a harmonic trap on preparation in a 
superlattice:- If a harmonic trap Vt rap is present for the 
duration of the preparation, we find that the character 
of the final state in terms of the pair momentum distri- 
bution is close to the 77-condensate, though the peaks are 
slightly broadened and the density profile will correspond 
to a trap. This state is close to an excited eigenstate in 
the presence of the trap, and for U 3> J is well approx- 
imated by an ansatz rf A — J^i Aic\ *c| , , where A e M M 
correspond to the ground state wavefunction of a single 
bound particle with tunnelling amplitude — J 2 /U in the 
presence of a trapping potential 2 Vtrap ■ 

Revealing and characterising imperfections with an 
77— condensate:- The 77-condensate is an exact excited 
eigenstate of the Fermi Hubbard model, and the cor- 
relation functions will be both sharply peaked and sta- 
tionary, unless there are errors in the state preparation 
or implementation of the Hamiltonian. Broadening and 
time-dependence of the correlation functions can be used 
to reveal imperfections, and also to characterise their 
source. We consider an initial eta-state, with delo- 
calised impurity atoms (see below for more details), and 
in Fig. 4 we plot the time dependence of the height of 
the peak in the pair momentum distribution. In Fig. 4a 
we consider only the additional atoms, and in Fig. 4b we 
add also a weak additional harmonic trapping potential. 
As in Fig. 3d, increasing Ni reduces the height of the 
77-pairing peak. However, provided U > 4J, the result- 
ing pairs, and the correlation functions are stable as a 
function of time. For U < 4 J (not shown) the pairs can 
decay through collision with unpaired atoms 0], and 
the peak in the pair momentum distribution also decays. 
On the other hand, additional potentials will dephase 
the state, and cause decay of the peak, as shown for a 
very weak harmonic trapping potential in Fig. 4b. The 
rate of decay is larger for stronger traps due to faster 
dephasing, and unlike the effect of missing atoms, is in- 



dependent of U/J. This difference could be used in an 
experiment to characterise the source of defects in the 
final state. Note that in order to make this discussion 
independent of the form of the ramp, we have obtained 
the results in Fig. 4 beginning from a state of the form 

\l,N,Ni) = E W){J} = nti(-l)^Ln£ie*Ha : 
|0) where : . . . : denotes the ordering operator by site, i.e 

' y I := VxCy 1 if x < V an d — c y ivl otherwise. Note 
that \h,N,Ni = 0} = |?77v)- 

Outlook:- The preparation of the 77-condensate offers a 
testbed to verify the emulation of many-body Hamiltoni- 
ans in optical lattices, providing both a sensitive means to 
validate the implementation of the Hamiltonian, and also 
an important test case for state preparation schemes in- 
volving adiabatic ramps. These schemes are particularly 
important in light of the current experimental challenge 
to reduce entropies in order to generate states such as 
an anti-ferromagnetic phase of the Fermi-Hubbard model 
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